
PRELIMINARY SHUTTLE STRUCTURAL DYNAMICS 
MODELING DESIGN STUDY 


W 


GRUMMAN 


^ 73-/3 7 // 

NASA CR 112196 


PRELIMINARY SHUTTLE STRUCTURAL DYNAMICS 
MODELING DESIGN STUDY 


Prepared Under Contract NAS1-10635-4 by 
Grumman Aerospace Corporation 
Bethpage , New York 11714 

NOVEMBER 1972 


for 

Langley Research Center 
Hampton, Virginia 23365 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


FOREWORD 

The work presented in this Report was the combined effort of 
TRW Systems Group Incorporated acid Grumman Aerospace Corporation. 
Sections 2, 3, and 8 were prepared by TRW (A. E. Galef, Dr. P. B. 
Grote, and D. H. Mitchell respectively). Sections 4 and 6 were 
prepared by J. B. Smedfjeld, and Section II was prepaired by R. L. 
Goldstein following suggestions by S. Goldenberg. The remainder 
of the report and all the revisions and addition of updated 
information were prepared by M. Bernstein. 

The design study of the l/lOth scale model details was 
conducted by A. P. La Valle , F. L. Halfen, and P. J. Coschignano. 

This work was administered by the Dynamic Loads Branch, 

Loads Division, NASA Langley Research Center, Hampton, Virginia. 


This report has been approved by: 



Program Manager 



TABLE OF CONTENTS 


Section Page 

Introduction and Chronology ......... 1 

1 Rationale for Model and Proposed Model Features:. . . . . . . 1-1 

2 Requirements for POGO Analysis 2-1 

3 Requirements for Control System Design .... ,3-1 

4 Ground Wind Loads Analysis * . . . .. . . .... 4-1 

5 Lift-Off Loads Analysis . 5-1 

6 Gust-Loads Analysis ........ 6-1 

7 Engine Induced Response ...... . . 7-1 

8 Staging ... 6-1 

9 Test Requirements .... . .... . . , . . . . . . . , . . . 9-1 

10 Summary of Requirements Study . . . . 10-1 

11 Modal Coupling Applied to 1/15 Scale Model Test Data . . . . 11-1 

12 Design of 1/10 Scale Preliminary Model of Series -Burn 

Configuration 12-1 

References . ..... , . . . . . R-l 

Appendix A-Gust Response Analysis . . . . . . . . . . , . . . A-l 

Appendix B-l/10 Scale Dynamic Model Preliminary Stress 

Report . . . . . ... , . . . . . . . . B-l 


INTRODUCTION AND CHRONOLOGY 


INTRODUCTION AND CHRONOLOGY 

Structural dynamics problems during the Apollo program were very important 
and required considerable technology development concurrent with the design 
and flights of the various vehicles. Problems caused by POGO in particular 
required a high level of effort in very short time span. Increased analytical 
and test capability has since been developed and should help to reduce the 
impact of these problems on the Shuttle program. However, the Shuttle 
configuration will have .more complex structural dynamic characteristics than 
previous launch vehicles primarily because of the high modal density at low 
frequencies and the high degree of coupling between lateral and longitudinal 
motions. A preliminary structural dynamics model will provide for the early 
study of many structural dynamics problems, a means of evaluating the accuracy 
of the structural and hydroelastic analysis methods on a shuttle-like test 
vehicle, and a means for efficiently evaluating potential cost savings in 
structural dynamic testing techniques. 

This task was originally established in June 1971 to determine the optimum 
configuration for such a preliminary model. The original requirements for the 
study included the following: 

(1) Determination of requirements based upon a review of significant 
dynamics analyses (POGO, Control System Stability, Dynamic Loads). 

( 2 ) Selection of a generalized model configuration 

(3) Study of two candidate scale factors 

(4) Study of initial structural details and material substitutions 

(5) Preparation of Test Plans and Schedules 

The requirements study based upon the Shuttle configurations under 
consideration at that time was completed and reviewed on August 12 , 1971 . 

Copies of the report were submitted to the Technical Monitor on September 20, 1971. 

As a result of changes in the Shuttle Configuration and at the request of 
the Technical Monitor, the effort on the remaining portions of the study; was 
suspended until more definitive "configuration 'information became available . . 


An additional study, authorized by Modification No. 1, to further 
evaluate modal component synthesis techniques was completed in November 1971, 
and submitted to the NASA Technical Monitor on December 13, 1971. This study 
verified that the analytical treatment of measured data was feasible for the 
existing 1/15 scale NASA Langly Shuttle Model, and that further analysis was 
warranted. As part of the same effort, a review of the measurements required 
to compute a damping matrix from test data was conducted at TRW. The procedure 
proposed is described in an AIAA Journal note "Method for Constructing a Pull 
Modal Damping Matrix From Experimental Measurements" by T. K. Hasselman, 

April 1972. The attempt at using this method on available NASA and TRW test 
data showed that the presence of off -re sonant modes in the responses made the 
data unsuitable, and additional processing was required. This effort was 
transferred to Task NAS 1-10635-8 for processing vibration test data. 

Modification No. 3 authorized a change in the basic objectives of the 
task by directing that a general arrangement and some structural details of a 
1/10 scale model design for the series-burn pressure-fed booster, which had 
been selected as the baseline configuration, be developed. In accordance with 
this direction, model structural studies and an accompanying stress analysis 
was undertaken starting in January 1971. Proposed details of the LO^ tank, 
intertank skirt, and LH^ tank were discussed with the NASA Technical Monitor 
in January 1971. The preliminary design and analysis was completed and 
submitted on March 31, 1972. 

This completed the model design effort under this Task Order. The' design 
details for the tanks and orbiter developed during this effort are presently 
being used as part of the design of a complete model under NASA 1-10635-11. 

An additional brief study of the suitability of the proposed model design 
details for obtaining structural dynamics data for use in verifying POGO 
analytical methods was completed by TRW in- August 1972. . 

All the work accomplished under this Task Order, as described above, has 
been included in this report. Several of the sections have been revised 
slightly from their original submission to include additional available 
information. 
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1. RATIONALE FOR MODEL AND PROPOSED MODEL FEATURES 
1.1 Necessity for a Preliminary Shuttle Dynamic Model 

Structural dynamic effects will be much more important for the Shuttle 
than for past launch vehicles because the configurations evolving from the 
various studies all possess large numbers of low-frequency modes with con- 
siderable amounts of coupling between motions in orthogonal directions. These 
effects are analytically accounted for by using mathematical models derived 
from a large-scale finite -element representation. Experience on less complex 
vehicles such as Saturn, LM, and F-lU has shown instances where important 
differences occurred between analysis and test results, for example, between 
the calculated and measured mode shapes and frequencies. It is imperative, 
therefore, that shuttle analytical procedures be checked (and modified or 
supplemented as required) by comparing typical results from these procedures 
with corresponding test data from a structure containing as many of the 
significant prototype features as is feasible. A preliminary dynamic model 
provides the only means of obtaining this complementary experimental data on 
a timely and cost-effective basis. Specific contributions of the dynamic-model 
program are: 

• Enhancement and verification of analytical modeling procedures in areas 
where these are reasonably well developed (e.g., definition of finite- 
element structures models and lumped mass dynamics models) 

• Complementary empirical data in areas where analyses are relatively 
weak (e.g., representation of fluid/tank-wall interactions for highly 
flexible thin shell tanks, or modeling of partially buckled structural 
elements) 

® Investigation of effects of system nonlinearities (e.g., preloaded 
interstage restraint systems) 

In addition to these analysis-related, requirements for a dynamic model, 
some further benefits include: 

• Early verification of important preliminary vehicle dynamic character- 
istics of a typical shuttle design 

• Disclosure of previously unanticipated dynamics problem areas 
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• Definition of optimum configuration changes for eliminating critical 
problem areas 

• Development and demonstration of efficient full-scale prototype testing 
procedures. 

In summary, while fulfilling its primary mission of providing vitally- 
needed supporting data in key analysis areas, the model will also furnish 
means for investigating critical design and testing problems. 

1.2 Rationale for Selecting Configuration Features 

The primary purpose of the preliminary dynamic model is to verify the 
procedures used in analyses. Therefore the major emphasis in this limited 
study establishing model requirements was to review each of important analytical 
procedures in detail and determine the necessary structural dynamics information. 
In most instances this consisted of determining the modes in which major input 
forces and responses occurred, and the significant structural details in these 
modes. The frequency ranges required were determined principally from analyses 
conducted during shuttle configuration studies or from Saturn V results. 

j 

The model is also required to incorporate as many significant typical 
prototype features as possible so as to reveal the maximum number of problem 
areas. In order to limit costs, however, these features should be simplified 
and idealized while still simulating the dynamic behavior. 

1.3 Summary of Recommended Configuration Features 

The significant structural features required for the major elements of a 
model to be representative of a prototype for each analysis reviewed and most 
likely to reveal problem areas are summarized in Section 10. Based upon these 
considerations, it is recommended that a preliminary dynamic model incorporate 
the features listed below. - 

• Orbiter 

The orbiter fuselage should be represented by three sections consisting 
of a forward portion containing the. crew area and forward payload . ... 
attachments, a mid-section containing the wing carry-through and the 
aft payload attachments, and an aft section containing the engine and 
fin support structures. The crew compartment should be a simple closed 
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shell containing extensions of the two side longerons with provisions 
for attachment of weights representing equipment. The mid section 
should have removable payload doors which have torsion carrying 
capability but add no bending stiffness. The basic U shape frame for 
the mid section should be used along the entire fuselage. 

• External Tanks ( 

The external tank should consist of a representation of an LO^ tank 
designed to hold water, an intertank skirt with provisions for attaching 
SRM's, an LH^ tank with internal frames and external fittings to attach 
the Orbiter, and an aft skirt with provisions for attaching the booster. 
The LO^ tank should be represented by a monocoque shell of uniform 
thickness.. The bottom dome should be of uniform thickness with attach- 
ment for measuring pressure and for drain and fill, capability. No 
provision is required for slosh baffles or simulated feedlines. 

Internal pressurization can be limited to a low value (2 psi). 

The intertank skirt should be of monocoque construction with no chem 
milling. Any frames or struts used to distribute loads internally 
should be represented by the proper geometry, however, in order to keep 
the design simple, the members should be designed to avoid buckling 
without requiring added stiffening. Interstage attachment members may 
be of simple cross-sectional shape rather than scaled down prototype 
values. Fittings may be oversized and simplified to limit costs. 

The LH^ tank should be of monocoque' construction with a minimum of 
chem milling. Skin gages should be heavy enough to avoid buckling; 

' t 

under model handling loads without requiring any additional stiffeners.'. 
The internal frames which transmit the orbiter loads should have the 
proper geometry but the members may be of simpler cross sectional shape 
and should be designed to avoid buckling without stiffening. Fittings 
may be oversized and simplified, to limit costs. 

The aft skirt should be of simple monocoque construction, similar to 
the inter-tank skirt described above. 
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POGO CONSIDERATIONS IN SPACE-SHUTTLE STRUCTURAL DYNAMIC MODEL 
2 . 1 Brief Description of Analysis Procedures 

POGO is a self-excited oscillation that results from coupling of 
structure, propellant, the pumps and the motors of a rocket propelled vehicle. 

Qualitatively, the POGO phenomenon occurs because vehicle accelerations 
and associated accelerations of the propellants in tanks and feedlines cause 
pressure changes in the propellant at the pump inlet. These pressure changes 
result in changes of net thrust of the rocket motors. The phase of these 
thrust changes is usually such as to cause structural accelerations in the same 
direction as those initially postulated. There is then "positive feedback," 
and the system is potentially capable of self-excited vibration. 

A determination for a particular configuration of a specific vehicle 
of whether POGO is likely to occur, may be made from an evaluation of the 
closed-loop eigenvalues of the coupled engine -structure -line system or, 
alternately, by a calculation of the open-loop frequency response of the same 
system with the loop opened at pertinent points. Regardless of whether open- 
loop or closed-loop analysis is selected, it is necessary to enter the same 
parameters in the analysis. These parameters are: 

(1) Structural response, as functions of magnitude and frequency of 

exciting force at engine locations 

a. Engine motions 

b. Pump motions 

c. Motions of fluids in fuel and oxidizer tanks, or alternately 

d. Pressure of fluids at tank-feedline interfaces 

(2) Feedline Characteristics 

a. Pressure at pump inlets caused by pressure at tank-feedline 
interfaces and motion of lines and pumps*- 


*This pressure is often influenced strongly by the naturally occurring or 
artifically provided compliance at the pump inlet, and is therefore determined 
separately from Item 1-d, above. 
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(3) Pump-Engine Characteristics 

a. Variations in net thrust caused by variations in pressure and 
flow rates at the pump inlet 

b. Compliance at pump inlets 

All these parameters may be determined, with more or less accuracy and confidence, 
by analysis at early stages of the vehicle design. But they must all be test- 
verified prior to flight in order to improve confidence in results of POGO 
stability calculations. 

The structural characteristics should be test-verified as early as possible 
because of the tremendous schedule and cost impact that could result from 
unforeseen structural characteristics that are sufficiently unfavorable to 
indicate a need for significant changes. A structural dynamic scale model, 
which may be considered as a hybrid between analysis and a full scale structure, 
can be of great help in confirming the analytic techniques used to calculate 
structural characteristics. 

The feedline characteristics, as functions of the compliance prevailing at 
the pump inlet, have been considered in the past to be the characteristics 
which are most amenable to analysis and therefore to require the least test 
verification. However, the convolutions of presently proposed space-shuttle 
feedlines are such that confidence in analysis is reduced, and test-verification 
prior to building of actual hardware will be appropriate. Models of feedlines 
will be valuable, but these need not be incorporated in the overall model for 
the reasons discussed later in Section 2.4.3.' Instead, the effect of the lines, 
when known from analysis and test can be incorporated in the overall mathematical 
model by reliable modal synthesis techniques. 

The pump and engine characteristics, including the cavitation compliance 
at the pump inlet which has been stipulated often to have strong effects upon 
the feedline characteristics, are dependent upon both fluid velocities and 
accelerations, which cannot be scaled simultaneously. Therefore, these 
characteristics cannot be derived from scale models . and_a scale .model cannot 
be used to provide a complete determination of the POGO stability of the space- 
shuttle. But because a structural dynamic scale model can provide some 
important information in a timely manner, it can. be. a very valuable tool for 


determination of POGO stability. Such a model is recommended. 

In the remainder of this section, considerations entering into the cost- 
effective and timely fabrication of such a model are discussed. Specific 
recommendations are made when possible, and areas requiring further studies 
are identified. 

2.2 Tanks and Fluids 

2.2.1 Tank Structure 

Tank bulging influences all longitudinal modes, and has a direct effect 
on the modal pressures at the tank-feedline interface. 

Tank bulging is affected primarily by the density of the fluid in the 
tank and by the stiffness of the tank walls and bottom. The density of the 
walls and bottom is of negligible concern, as is the stiffness (bulk modulus) 
of the fluids. 

The stiffness of concern for thin-walled tanks as will exist in space- 
shuttle is primarily extensional stiffness; shell bending stiffness has purely 
local effects which are overwhelmed by the extensional stiffness. As a result, 
it will be permissible to model the tank walls and bottom with any material of 
appropriate thickness having stiffness (EA) equivalent to the composite 
stiffness of the proposed tank design. This permits the use of reinforced 
plastic for the model tank wall, and this may be particularly convenient if the 
full scale design should incorporate different circumferential and longitudinal 
stiffness. . An appropriate selection of amount of reinforcing material in the 
warp and woof of the . selected reinforcing fabric will permit adequate simulation 
of tank walls. , 

2.2.2 Fluids 

Present planning involves fueling the external tanks with liquid hydrogen 
and liquid oxygen. 

A mixture of water and alcohol could be made to provide precise density 
-simulation of liquid -oxygen j as -could many other combinations of miscible fluids'. 
Any such mixture would have substantially greater acoustic velocity than the 
highly compressible liquid oxygen but, as was noted,- this difference is not 
significant for tank bulging modes', 



The very low density of liquid hydrogen, which makes it difficult to 
simulate with a room-temperature fluid, usually makes it precise simulation 
unnecessary. Thus, when testing the "Apollo" full scale Dynamic Test Vehicle 
it was conclused that negligible accuracy would be lost if the liquid hydrogen 
were replaced in the scale model by an equivalent weight of water. An equivalent 
conclusion may well be reached for the proposed space-shuttle model. 

If the lack of need for simulation of the liquid hydrogen cannot be 
demonstrated for space -shuttle, then it will be possible to provide reasonable 
density simulation with a uniform grid of closed-cell plastic foam blocks, 
constrained against floating to the surface, displacing appropriate amounts of 
water. Appropriate density plastic beads, either foam or phenolic "micro- 
ballons," might prove to be acceptable without any fluid. 

Finally it is noted that while simulation of fluids by spring masses will 
be useful in some model studies, particularly those involving slosh modes where 
gravity scaling is impossible, such simulations do not appear to be either 
necessary or sufficient for a model designed to yield POGO information. 

2.3 Structural Transfer Functions Required 

The structural characteristics which have been hoted as' necessary for 
quantitative POGO analysis have been stipulated to be frequency-dependent 
functions. As a result, the structure is most conveniently described for POGO 
analysis by sets of characteristics of its normal modes. The pertinent 
characteristics are: 

(1) Modal frequency 

(2) Modal damping . 

(3) Generalized modal mass • 

(4) Modal motions of engines, pumps, feedlines, and propellant in tanks, 
normalized to the same point used for normalization in calculation of 
the modal mass 
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The modes of concern are those having substantial engine motion in the direction 
of engine thrust, and having modal frequencies near oxidizer* feedline resonant 
frequencies. In the case of the External Tank and SRM configuration, the 
frequency of the lowest feedline resonance of concern may be as low as 2.5 Hz, 
and only those structural modes with frequency less than approximately 25 Hz 
need be considered. For the Orbiter and External Tank configuration, or for a 
proposed configuration which has oxidizer tanks aft, however, it will not be 
possible to define a usefully low "upper limit of concern" until after accumula- 
tors have been assigned to the pump inlets, and reasonably sized. The "upper 
limit" for the orbiter, or for a booster with oxidizer tank aft, may be higher 
than 30 Hz. 

2.3.1 Structural Details 

Low "maximum" frequencies in the actual vehicle are desirable from the 
viewpoint of designing and building a model, for only those local assemblies 
having subsystem frequencies less than approximately one octave above the 
"maximum" frequency need be modeled in detail. Stiffen subassemblies are 
adequately modeled when their rigid body inertial properties are reproduced in 
scale. Conversely, some structure will be adequately modeled if its stiffness 
alone is reproduced in scale; an example would be the interconnecting structure 
between the booster and the orbiter. 

2. 3.1.1 Thrust Structure 

The thrust structure, which interconnects the engines and attaches to the 
orbiter, may be a dominant factor in POGC- significant modes if the maximum 
frequency of concern extends reasonably high, and plans should be made for 
careful modeling of that structure. The expected simplicity of the orbiter 
thrust structure is fortunate in light of the possibly substantially higher 
frequency regime of concern which may require valid simulation of modes as high 
as 30 Hz. This implied the. need for a scaled replica of orbiter thrust 
structure if the objective is a model with good fidelity to the prototype. 

*Fuel feedlines are. presently .expected to have little influence on POGO sensi- 
tivity of the space shuttle because of the relatively shorh lengths of the 
fuel feedlines and because of the mass flow rate of liquid hydrogen will be of 
the order of 20$ the flow rate of the liquid oxygen. These factors make large 
reductions in thrust oscillations for fuel side velocity changes as against 
oxidizer-side oscillations, and therefore make fuel natural frequencies of less 
concern. 
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Unless oxygen feedline accumulators should be firmly incorporated in the design, 
it is recommended that plans be made for including a scaled replica of this 
subsystem. 

2. 3. 1.2 Lines 

An important subassembly, whose modal frequency is by definition within 
the frequency regime of concern, is the oxidizer feedline system. This is 
discussed separately in Section 2.4, following. 

2. 3. 1.3 Joints 

Good design of joints usually aims for 100$- joint efficiency, which means 
that the joint will fail under the same load as skin adjacent to the joint. In 
order to approach the goal of 100% efficiency, it is generally necessary to 
weld a reinforcing ring to the skin, so that the local stiffness and weight of 
the joint will differ from that of adjacent continuous structure even when the 
structural efficiency reaches the desired 100 percent. 

Purely local increases in stiffness, even when quite large, have negligible 
effect upon structural modes. If analysis shows that any of the doublers 
hypothesized above increase local stiffness, their weight effects may be accounted 
for by adding some local weight to the model, and their stiffness effects may 
be neglected. 

Local reductions in stiffness have more effect on modal characteristics 
than do stiffness increases, but the effects remain small in the absolute sense 
unless the stiffness reduction is large and; this is not likely to occur in the 
case of a normally designed joint. If, however, calculations indicate that the 
stiffness reduction is indeed large, this can be simulated by cutting away a 
small region of skin much more easily and economically than by attempting to 
reproduce the joint. 

2 .4 Propulsion Study 

Transfer Function Requirements 

The, transfer- functions -which -are -reqiiiredrwith' '.the structural transfer 

functions to complete the POGO loop are: 
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'• 3V 3P sO> 3 V P sf> 3P c /3i, DO’ 3 VV 

Z ' . Vt’V.-' Vt ■ 

Where P c is engine chamber pressure 

P , P , are inlet pressures of oxidizer, fuel, pumps, respectively 

SL) S'X 

P^ is pressure at tank bottoms 

Wpf, is fluctuating flow rate in fuel, oxidizer, pumps respectively 
Y , Y , is motion of pump and line support points, respectively 

6 ii 

2.4.1 Engine 

The first set of functions, providing the chamber pressure and hence the 
thrust, are characteristics of the pumps and engines. As has been previously 
stated, these cannot be derived from scale . models. 

2.4.2 Lines 

The second set of functions, yielding the pressures at the pump inlets, 
are functions of the geometry and stiffness of the feedline structure, the 
density and bulk modulus of the propellants, and the location and stiffness of 
naturally occurring or artifically provided lumped compliance at any point of 
the line. Direct replica simulation of this set of parameters in model scale 
will be exceedingly difficult if possible , because: 

1. Line walls would have to have thickness as low as .002 inches even if 
the model were made as large as 1/10 scale 

2. The bulk modulus of cryogenic fluids tends to be very low in comparison 
to that of room temperature fluids 

3. Reynolds number (viscosity) considerations will probably enter, because 
of both the diameter reduction of a scale model and the necessary 
absence of "steady state" flow rate in a model 

Some of these considerations can be overcome by, for example, using a 
highly compliant material for pipe walls and thereby arriving at a reasonably 

A . ' 

correct overall "effective" bulk modulus, B,_ where ... - : 
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where B is the bulk modulus of the fluid 

D is local line diameter 
E is line wall stiffness, 
t is line wass thickness 

However, we recommend against such procedures and therefore against including 
simulated lines in the structural dynamic scale model because of the complexities 
and uncertainties of doing so. 

Additionally-, in deriving a mathematical model of a structure for POGO 
purposes, it customary to exclude the feedlines and their fluid from the basic 
structural model and later to couple the line modes with structural modes in 
the stability analysis. This has the important advantage of avoiding the need 
for calculation of new sets of structural modes for each change in feedline 
characteristics. 

Exclusion of lines from a structural dynamic scale model would have an 
equivalent, and cost effectwise, probably more cogent advantage since it would 
not be necessary to retest the overall scale model each time a different line 
configuration is considered. Based upon past experience in Titan and Apollo 
POGO studies, where improvement of POGO characteristics involved primarily 
changes to local lumped line compliances, this will be a great advantage indeed 
in the Space Shuttle model. 

The argument against including feedlines in the structural model is 
reinforced by experience on the 1/10 scale replica model of . Apollo-Saturn V. 

Lines were carefully reproduced to scale on that model, but no useful data was 
achieved from this inclusion because the model test procedures could not 
include any provisions to reproduce the effect of fluid flow and cavitation at 
the pump entrance. Consequently line modes did not interact with structural 
modes in a realistic manner and useful correlation between feedline pressures 
and forces and motions could not be attained.-. 

It will, however, be useful to build models of the lines alone, or a cold- 
flow facility which simulates the pump dynamics as well. . Such models- could 
help refine techniques for analytic calculations of line nodes, as functions of 
cavitation and/or Accumulator compliance, of the complicated line configuration 
presently planned. 
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A model of lines -alone would not be constrained by the scaling factor 
selected for the overall model, or by the material of the walls. Further, 
since the model would be for the purpose of verifying analytic techniques 
rather than achieving specific directly applicable numerical results, detailed 

I 

geometry would not. have to be reproduced. Thus a subassembly model with lines 
of sufficient diameter to facilitate making pressure measurements and to 
minimize viscosity effects could be built without the need to have straight 
sections of lines as long as the scale factor or diameters would require if 
detailed geometry were to be maintained. This would substantially facilitate 
testing. 

2.5 POGO Related Test Objectives and Data for Proposed Model Tests ^ 

2.5.1 Introduction 

The structural dynamic complexity of the Space Shuttle configuration 
poses potential problems of a greater magnitude than have heretofore been 
experienced on spacecraft/launch systems. To enhance the understanding and 
assess the problems envisioned, a structural dynamic scale model has been 
proposed upon which a variety of tests could be performed. In a very general 
sense the objectives of this model test program can be grouped in three broad 
categories as follows: 

• Enhancement and verification of analytic modelling techniques and 
procedures 

• Disclosure and identification of unanticipated structural dynamic 
interactions and/or coupling problems 

• Development and demonstration of appropriate testing techniques/ 
procedures and data acquisition/evaluation techniques applicable to 
full scale 'prototype testing. 

One of the prime outputs of the model test program will be the determina- 
tion of the modal characteristics of the scale model against which analytic 
models will be evaluated. The complexity of the Space Shuttle configuration 

(l) This section was added in August 1972 after, the proposed model configura- 
tion for the parallel -burn Shuttle (with SRM's) was reviewed. 
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may require that the analytic modal characteristics be determined using 
component mode techniques, wherein subsystems such as the solid rocket motors, 
the H/O tankage, and the orbiter are individually and separately analytically 
modelled and then combined to generate shuttle system modal characteristics. 

An overall model test/analysis program as generally depicted in Figure 2-1 is 
envisioned. In support of such a program, the proposed design of the scale 
model H/O tankage (Reference 2-1) has been examined relative to the determination 
of modal tank bottom pressure coefficients necessary for a POGO analysis of the 
Space Shuttle . 

2.5.2 Background 

In longitudinal stability (POGO) analyses, an important set of parameters 
can be the "tank bottom pressure coefficients" (Reference 2-2). Such 
coefficients relate the oscillatory fluid pressure in the tank at the outlet 
to the feed system to a reference modal displacement or acceleration. The 
resulting tank bottom pressures are applied as a driving force to the feedline 
providing one of the coupling terms between the structural and propulsion 
systems in the closed loop interaction (POGO). The proper determination of 
these terms has only recently come under investigation as launch vehicle systems 
have become more complex (References 2-3, 2-4, 2-5). 

Reference 2-4 gives a summary of dynamic modelling technology for liquid 
propellant launch vehicles circa 1969-70. At that time, approximate methods 
were developed to more adequately represent the interaction between liquid and 
structure in the propellant tanks. These techniques were sufficient for 
calculating the longitudinal modes of the vehicle, but Reference 2-4 states, 

"The method is not applicable where higher tank modes or local liquid effects, 
such as tank bottom pressures, must be predicted accurately." 

Accordingly, an improved technique was used (Reference 2-3) to calculate 
tank bottom pressures for POGO analyses. The state of the fluid continuum is 
defined in terms of a velocity potential and finite elements are used to model 
the flexible tank. Surface wave effects are neglected and the liquid surface 
is assumed to remain flat. The analysis was compared with model tests but the 
experimental data were rather limited. ' 
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Another finite element technique (Reference 2-8) is similar to that of 
Reference 2-3 in that the fluid is treated as a continuum, but the problem is 
formulated in terms of an integral equation and discrete elements are used on 
the fluid boundaries. By using washer-like, annualr elements on the free 
surface, axisymmetric liquid sloshing effects are included (i.e., surface waves). 

After the Saturn V experienced POGO on the S-II stage, vibration tests 
were conducted on an S-II Lox tank (termed the mini-stage). These experiments 
were compared with an analysis to determine tank modes and tank bottom pressure 
coefficients (Reference 2-5) which used the above mentioned finite element 
program. This recent paper (Reference 2-5, January 1972) presents a comparison- 
between test and analysis. Although there are only limited tank bottom pressure 
data the modes and pressures compare reasonably well with the computer model. 

In prior scale model testing of both the Titan III (Reference 2-6) and 
Saturn (Reference 2-7), the need for tank bottom pressuee coefficients was not 
compelling and no significant information was obtained although an attempt was 
apparently made to determine tank pressures. A similar situation prevailed 
during the full scale testing of Saturn V and no useful pressure information 
was. obtained from these tests (Reference 2-4). 

Presuming that the Shuttle scale model test program will contain, in some 
sense, the elements as depicted in Figure 2-1, the following sections provide 
objectives and data requirements pertinent to H/O tankage tank bottom pressures 
for POGO analyses. Similar examinations should be made for the other subsystems 
(solid rocket motors, and orbiter) and for the other dynamic analyses (controls, 
transient loads, separation dynamics, etc.). 

2 . 5.3 Test Objectives • 

These objectives pertain to the H/O tankage as a component and as a part' 
of the total Space Shuttle system, and are oriented specifically toward obtaining 
data necessary for POGO analyses. 

For the H/0 tankage alone, as a Space Shuttle subsystem or component, the 
prime test objective -is: ^ 


(1) The determination of L0 p tank bulging/aft dome longitudinal mode 

shapes and frequencies and corresponding modal, tank bottom pressures 
under several propellant loading conditions. Modes up to approximately 
20 Hz (prototype) with the tankage excited at the LO^ tank "Y" ring 
are required. 

For the H/0 tankage integrated with the other subsystems (orbiter; with 
and without solid rocket motors) the objectives are: 

(1) Determine the free-free coupled modal characteristics of the total 

Shuttle system (up to 20 Hz prototype) under several propellant 
loading conditions. - 

(2) Determine tank bottom pressure coefficients for each of the system 
modes identified. Both LO^ and LH^ tank bottom pressure coefficients 
should be determined, but any compromises should be resolved in favor 
of the LO p tank measurements. 

2.5.4. Data Requirements 

For the H/0 tankage as a component or subsystem both accelerometer and 
pressure data are required. Placement of three rays of 4 accelerometers on 
the LO2 tank aft dome is a possible scheme to determine the first few dome 
mode shapes (Figure 2-2). The three rays of accelerometers will permit 
averaging around the dome to help define the modes and the lateral acceler- 
ometers above the "Y" will permit detection of tank bulging.' The intent is 
to define a dome mode shape that can be correlated with analytic results and 
to obtain a tank bottom pressure coefficient (psi/g or psi/in) for correlation 
with analytic predictions. 

For the Space Shuttle system model, both tank bottom pressure transducers 
along with sufficient accelerometers such that system modes up to 10-15 Hz - ' 
(prototype) can be defined are required. In particular modal data are required 
at the following locations for POGO analyses: 

• 3 axis accelerometers at the orbiter engine location(s). 

• 3 axis accelerometers at the L0 p feedline attach points to the 
structure (in both orbiter and H/o tankage) 


• both LCL, and LH 2 tank bottom pressures and accelerations 

• 3 axis accelerometers at the turbopump location(s) (if different from 
the engines) 

Pressure transducers and accelerometers should have flat phase and 
amplitude response characteristics in the range of approximately 20-300 Hz. 
Oscillatory pressures will probably be less than 5 psi. (Static pressure in 
an unpressurized LO^ tank full of water is ~ 2.5 psi.) The detail of the tank 
bottom deflections may not be required if good correlation can be obtained in 
component or subsystem tests. 

These data are required at several propellant loading configurations a 
minimum of which are (assuming an approximate 350 second powered flight and a 
100 second solid rocket motor burn) : 

t = 0 Lift off 

t 100 Booster burnout (booster attached) 

t =« 100 Booster burnout (booster released) 

t 350 Main engine burnout 

All data should include amplitude, phase and frequency, as a function of input 
amplitude and frequency. 

2.5*5 Configuration Considerations 

Since one of the prime objectives is the analytic/test correlation 
relating tank bottom pressures and modal deflections (or accelerations) rather 
than obtaining directly scaleable results to the prototype, several cost 
effective compromises seem appropriate. Tapered thickness of the lower bulk- 
heads is not necessary for the test/analytic comparison as long as the bulkhead 
modes remain in the proper frequency range relative to the overall system modes. 
Local stiffening and mass effects such as the tank outlet fittings/line or 
reinforcement structure should be included in the model. Water can be used as 
an L0 2 simulant. Light weight plastic spheres have been used before to simulate 
the mass distribution of the LH 2 (Reference 7), but the adequacy of this 
simulation of a liquid for obtaining tank bottom pressure measurements is 
unknown. Because the LHg propellant and feedline dynamics are not expected to 
contribute significantly to potential POGO problems due to the low density of 
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LH 2 , the uncertainties resulting from the use of such spheres will probably 
not be important. However, this choice is probably adequate for the overall 
system modal definition.- 
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Figure 2-2 Tank Bottom Instrumentation 



CONTROL SYSTEM CONSIDERATIONS 


3. COMROL SYSTEM ANALYSIS 

3.1 Brief Description of Analysis Procedures - 

The shuttle flight control system design necessitates an investigation of 
the interaction among the shuttle rigid and flexible body dynamics, the 
control system sensors and actuators and the external disturbances to which the 
vehicle is subjected. The coupling between the : controls, dynamics and the 
disturbances is described analytically via a closed loop formulation of the 
equations of motion such as depicted in Figure 3.1. The stability of the 
control system and verification of the system performance in the t im e domain 
must be established for each 'of the many flight maneuver's of shuttle vehicle. 
These include lift-off, normal atmospheric and orbital flight, staging, and 
aboil maneuvers which may occur throughout the vehicle's flight path. The 
performance of these analyses will begin with simplified analytical models 
(i.ev, a rigid body depicting the vehicle dynamics subject to a linearized 
control system) and culminate with complex models which include the flexible 
degrees of freedom coupled into the rigid body motion and including the highly 
non-linear control system. 

The success of the advanced control system analysis of the shuttle vehicle 
depends in part on obtaining a realistic description of the vehicle's structural 
dynamics. The size and complexity of the proposed shuttle structure requires 
that much of this information be derived from sub-scale model testing. Some of 
important aspects of the controls -dynamics interaction of the shuttle vehicle 
are discussed individually below. , * . 

The structural feedback into the control system affects both sensors and 
actuators. Sensors (rate gyros and accelerometers) are primarily intended to 
sense rigid body rotations and accelerations, but their output includes errors 
due to structural bending rotations and linear flexible body accelerations. 

These flexible body effects are minimized by judiciously placing the sensors 
so that the local structural dynamics response will on the whole be minimal, 
and by using several sensors placed at various points on, the system such that 
a comparison or averaging of these sensor outputs by onboard computers allows . 
an accurate prediction of rigid body rotations and accelerations. 


In order to evaluate sensor error, accurate modal data (modes, frequencies, 
generalized masses) are required. Secondly, the existance of local modes must 
be examined to ensure that the gross vehicle modes provide a sufficient basis 
for error processing. In the same manner, the flexibility of the points of 
action of control actuators must be considered in the analysis of structural 
feedback. 

Control actuators include booster and orbiter engines which may be 
gimballed or have variable thrust, attitude control jets and aerodynamic 
control surfaces. Engine mount point and engine actuator flexibility must.be 
considered as well as the flexibility of -points of attachment of control 
surfaces. 

Propellant slosh effects on the control system are of major importance due 
to the predominance of propellant fluids in the total mass distribution of the 
shuttle vehicle. Slosh modes may interact with one or more of the control 
system frequencies during both steady state and transient maneuvers. . The 
slingshot effect is an example of the latter category. It may take place at 
the end of an engine burn when energy stored in the deformation of fuel tanks 
during high acceleration periods is suddenly released causing liquid residuals 
to be propelled towards the forward bulkhead at high velocity. 

Other mechanisms of structural feedback involve aeroelastic effects which 
may couple the rigid and flexible body dynamics and localized deformation due to 
the mass and inertia distribution of major components. 

A sub- scale model of the shuttle vehicle can be readily used to obtain the 
gross natural modes of vibration within the control system band-pass (~ 7 Hz). 
Yet, the complexity of struct lire makes it very difficult and costly to attempt 
a replica scaling of local flexibility and structural damping effects which, 
as was noted above, play a major part in the controls-dynamics interaction. 
Furthermore, some phenomena, such as fuel sloshing are subject to different 
scaling laws than the natural modes of vibration (i.e. , slosh frequency 
typically is scaled by the square root of the length ratio) and thus are 
difficult to scale accurately. 


Fortunately, replica scaling representing identical geometric scaling 
and using prototype material throughout, the model is not necessary to 
establish the validity of the analytical simulation of the dynamic interaction 
of the shuttle vehicle. Instead, general geometric scaling may be employed, 
where the localized flexibility and a number of slosh modes are expressed in 
terms of . "distortion parameters" that indicate the degree of distortion of 
the individual parameters used. An analytical simulation of the closed-loop 
controls-dynamics interaction of the sub-scale model is then performed and the 
results are compared -with a laboratory test performed on the scale model. Good 
agreement between the analytical and the test results permits extension of the 
simulation, using the aforementioned "distortion parameters" whenever possible, 
to cover the shuttle flight configuration. 

3.2 Structural Information Required 

The structural dynamics data required or validated from the sub-scale 
model of the shuttle vehicle are as follows: 

(1) Lift-Off Configuration 

(a) First three natural frequencies and modes of vibration of the 
cantilevered booster/orbiter combined structure 

(b) Modal displacement at engines and accelerometer locations 

(c) Modal rotation at the mount point of each rate-gyro 

(2) In-Flight Configuration 

(a) Natural frequencies and modes of vibration of the free-free '• 
booster/orbiter combination and the'orbiter alone within a 
frequency band up to 7 Hz 

(b) Modal displacement at engines and accelerometer mount point 

(c) Modal rotations of the mount point of each thruster and rate gyro 

(d) Localized stiffhess coefficients of major components measured 
statically 

(e) Structural damping coefficients of all major modes of vibration 

(f) - Discrete damping at constrained hinges 

(g) Slosh modes and frequencies up to 7 Hz 


These data should be obtained for several vehicle configurations including:. 

• Booster/orbiter configuration with booster propellant loading of 
100$, 75$, 50$, 25$, and 0$ 

• Orbiter configuration with propellant loading of 100$, 75$, and 50$, 
25$ and 0$ 
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Figure 31 . Block Diagram of' the Shuttle Control Loop 
with Flexible Body Dynamics 
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GROUND WIND LOADS CONSIDERATIONS 


4. GROUND WIND LOADS ANALYSIS ■ 

4.1 Brief Description of Analysis Procedures 

Dynamic loads due to ground winds may be caused by one or more of several 
phenomena: responses to turbulence -induced unsteady drag forces, oscillatory 

lateral forces due to vortex shedding, galloping instability, and stop-sign 
flutter. 

Turbulence -induced loads may be determined analyticailying using a model 
of the environment such as that contained in Reference (4-1) together with 
vehicle distributed drag data and structural dynamic characteristics. The 
approach used would be a cross -spectral analysis such as that described in 
Reference (4-2). Since a shuttle-type vehicle could be excited by turbulence 
in both translation (bending) and roll (torsion) , antisymmetric as well as 
symmetric modes would have to be considered. 

An alternate scheme for obtaining ground -wind gust loads would be to 
utilize elastic models in a wind tunnel which includes a simulation of the 
atmospheric turbulence as discussed in Reference (4-3). Power spectra of the 
shears, bending moments, and accelerations would then be obtained directly from 
the models. For more' detailed design information than might be available from 
the model instrumentation, the measured results could be compared with the 
analytical results obtained in accordance with the procedure of the previous 
paragraph; these comparisons at selected locations could then be used as factors 
to be applied to the remainder of the analytical results. : It should be noted 
that this subject is still in the developmental stage, and difficulties in 
reproducing the desired spectra and spatial correlations may be significant. 

4 . 2 Structural Transfer Functions Required 

Analyses and/or wind-tunnel tests to determine dynamic loads due. to 
ground winds will require a knowledge of the cantilever modes of the shuttle 
restrained on the launch pad by the hold-down structure. As noted above, both 
symmetric and antisymmetric modes are of interest. 
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Turbulence -induced loads will occur primarily at low frequencies - in 
part due to the spectral distribution of energy in the atmosphere, in part 
because the vehicle restraint prevents the alleviation of very-low-frequency 
gust inputs by rigid-body vehicle motions, and in part due to the partial 
cancellation of positive and negative generalized forces for mode shapes higher 
than the fundamental. Past experience with essentially cylindrical launch 
vehicles has shown that only the fundamental cantilever mode makes a significant 
contribution to the total mean- square response. In extending this experience 
to the space shuttle, the concept of a "fundamental mode" requires some 
generalization. First, as stated above, both symmetric and antisymmetric modes 
will have to be considered, i.e., bending normal to the wing plane, lateral 
bending, and torsion (roll). Also, since the wings will give rise to much of 
the aerodynamic loading, the fundamental wing modes, as -well as the fundamental 
cantilever fuselage modes may have to be considered. 

The phenomena of, vortex shedding, galloping instability, and stop-sign 
flutter all involve overall vehicle oscillations in the fundamental cantilever 
modes. Wing flexibility may influence the last two phenomena, but in any event 
the degree of dynamic similarity discussed in the previous paragraph on 
turbulence should also be adequate for these other considerations. 

The approach to determining loads due to vortex shedding will be largely 
experimental and/or empirical with several different variations being possible* 
The one which would utilize the most analysis would involve wind-tunnel testing 
of essentially rigid models to determine only the oscillatory lift forces. A 
power- spectral response analysis would then be used in conjunction with a 
structural dynamic mathematical model to determine dynamic loads. Another 
approach would be to test a dynamically similar wind-tunnel model to obtain 
responses and loads directly. Regardless of what type of model is used, 
however, careful attention must be paid to the effects of Reynolds number 
differences between the wind-tunnel and full-scale conditions. 

Since the shuttle configurations are symmetric only with respect to one 
plane, and since vortex .shedding could produce -lift -forces normal to “this 
plane, both symmetric and asymmetric modes will have to be considered. On the 
other hand, for the complex noncylindrical configurations that characterize the 
shuttle designs, the extent of vortex shedding is still in doubt. 
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The phenomena of galloping instability and stop- sign flutter are in one 
sense similar to vortex shedding and in another sense opposite to vortex 
shedding. All three phenomena must be studied primarily by experimental and/or 
emperical methods, but whereas the significance of vortex shedding is probably 
diminished in changing from a cylindrical to a shuttle configuration, the other 
two problems were not important previously and now may be very significant. 
Again, asymmetries in the vehicle configuration require that both symmetric 
and antisymmetric structural vibration modes be considered in any mathematical 
or physical model. 
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5. LIFT-OFF LOADS ANALYSIS 

5.1 Brief Description of Analysis Procedures 

Dynamic loaids at lift-off are calculated as the response of a vehicle, 
initially deflected by wind loads and engines, to the forces induced as the 
hold-down restraints are released. The initial deflected position can be 
determined from the response of the cantilivered vehicle to the combination of 
appropriate engine -induced forces, steady winds, and turbulence, obtaining an 
equivalent static deflection using the matching technique of reference (5-1). 

An elastic wind tunnel model can be used as an alternative means of determining 
responses to appropriate winds. 

For purposes of analyses the vehicle is represented by the free elastic 
modes j with initial conditions being described in terms of deflections in free 
modes plus rigid rotations. An alternate approach which in some cases is less 
difficult than matching deflections is to determine the reaction forces -while 
the vehicle is restrained by the hold-down and apply these as forces to the 
free modes. Transient forces representing the continued net difference between 
thrust build up and hold-down releases are then applied along with steady peak 
wind forces, and the maximum member loads are determined. The dynamic loads 
due to turbulence at lift-off are determined separately, with the flexible 
vehicle again represented by the free modes. The force applied in this case is 
the spectrum representing the horizontal wind turbulence. Again, the matching- 
condition technique of reference (5-2) can be used to convert member loads 
representing a 3-sigma response to the turbulence to a static load distribution. 
These results are then added to those due to lift-off transient and steady wind 
to determine total load magnitudes. 

5.2 Structural Characteristics Required 

The loads at lift-off are expected to designthe aft portion of the SRM. 
Dynamic loads at lift-off affecting members in this area would therefore be most 
significant. These would be determined by motions of the major mass items in 
the free modes. below 4 Hz, the most important probably being the liquid oxygen 
tank and solid rocket motors. The interstage flexibility and LOg tank to 
interstage flexibility are therefore important. 
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In this analysis the major forces are those applied principally in the 
axial direction through the vehicle engines and hold-down release system. The 
engine ignition start up transient is expected to have a relatively slow rise 
time. Current specifications limit the thrust buildup in each liquid rocket 
engine to a rate of less than 700,000 lbs/sec and the engines are expected to 
be started sequentially so that the time from 10 percent to full thrust should 
be relatively long. The solid rocket motor start-up time is also relatively 
slow, about .3 seconds. The hold-down release forces could be applied rapidly 
since the Saturn V system originally permitted release in .01 to .02 seconds. 
However, in the present system the release time is about 0.6 seconds, and a 
similar device might be used for the shuttle to limit dynamic response. The 
axial forces therefore will not contain many high frequency components, and 
the low frequency modes with appreciable axial motion at the engine thrust 
structure would be most significant . In the Grumman GUI configuration for 
example only 1 symmetric mode of the lowest 20 had appreciable deflection 
(above 0.1) at the SRM engine and hold-down location, and this occurred at 
2.5 hz. Many of the modes did have appreciable motion at the orbiter engines; 
however, the fluctuating thrust at higher frequencies will not add much to 
primary structure loads as mentioned in Section 7. 

Lateral forces are applied principally by the potential energy stored in 
the deflected cantilevered vehicle. The modes most significant in this 
context would be those in which static and dynamic response to ground winds and 
orbiter engines are highest. 

In order to obtain some insight into the frequency range corresponding to 
this criterion, the cantilevered modes at lift-off calculated for a 4.9 million 
pound parallel burn (orbiter/external tank/SRM's) configuration (Grumman 
Design Gill) were reviewed. A description of the mode shapes is listed in 
Table 5-1. The modes which are expected to be most significant are the 1st, 
3rd, and 4th, and 8th. The wing and fin modes for this vehicle are summarized 
in Table 5-2, and it may be noted that the fundamental frequencies are low and 
that the 2nd wing bending mode which is probably the upper limiting frequency 
for appreciable wind induced loads is at 8 hz. 
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CANTILEVERED SYMMETRIC MODES - GUI LIFT-OFF CONFIGURATION 



10 3.47 .934 M Fuselage Z-2nd Bend/ Z-2nd Cantilevered Z-lst Bend. Axial 

Nose-Z LIL, Axial Axial, X-lst Long. Pitch, Z-lst Bend 

Note: Rigid Modes denoted by Pitch, Axial, Plunge. Flexible by Z-lst Cantil.. T^no- 



TABLE 5-2 


. GUI ORBITER WING : AND FIN MODES 


FREQUENCY Hz 

RESTRAINT AT ROOT 

MODAL DEFORMATION PATTERN 

1.59 

Antisymmetric 

Wing 1st Bending 

1.8l 

Symmetric 

Wing 1st Bending 

3.5 

Flexible 

Fin 1st Bending 

3.92 

Antisymmetric 

Wing Tip (RCS Pod) Pitch 

4.45 

Antisymmetric 

Wing Tip (RCS Pod) Pitch 

4.96 

Antisymmetric 

Wing Tip (RCS Pod) Yaw 

5.01 

Rigid 

Fin 1st Bending 

5.73 

Symmetric 

Wing Tip (RCS Pod) Yaw 

. 7.53 

Antisymmetric 

Wing 2nd Bending 

8.27 

Symmetric 

Wing 2nd Bending 

9.52 

Antisymmetric 

Wing 1st Torsion 

10.18 

Symmetric 

Wing 1st Torsion 

11.51 

Flexible 

Fin 2nd Bending 

16.24 

Rigid 

Fin 2nd Bending 

18.31 

Flexible 

Fin 1st Torsion 

19.21 

Rigid 

Fin 1st Torsion 


Note ; Antisymmetric restraint allows no deflection at fuselage side 
wall and centerline. Symmetric restraint allows no deflection 
at fuselage side wall and zero slope at fuselage centerline. 
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6. GUST-LOADS ANALYSIS 


6.1 Brief Description of Analysis Procedure 

Space -shuttle flexible -vehicle gust -response calculations will generally 
be of three different types: (l) a power -spectral analysis for the maximum- 

dynamic -pres sure portions of the ascent and entry phases, of a mission using a 
gust spectrum representative of the high-frequency content of measured wind 
profiles; (2) a time-history ascent analysis in which the high-frequency content 
is retained in the wind profile; (3) a power spectral analysis for the near- 
horizontal-flight phases of the vehicle usage. Of the Various portions of a 
mission profile, the gust loads during ascent are expected to be the most 
severe. Also, during this flight phase, the structural -mode frequencies will 
be the lowest and the wing loading will be the highest. Thus, the contributions 
of the structural-mode responses relative to the rigid-body responses will be a 
maximum for this flight phase. Accordingly, in evaluating the importance of 
higher-frequency structural modes from a gust -response viewpoint, it was this 
flight regime that was considered. Details of this work are contained in 
Section 6.3, but the conclusion is that about 10 to 15 of the first 30 symmetric 
modes will have responses of sufficient magnitude in a gust-loads analysis of 
reasonable engineering accuracy. 

6.2 Structural Transfer Functions Required 

The structural modeling, from a gust-response viewpoint, should give 
emphasis to those modes which: (1) are readily excited by a gust normal to the 

vehicle, i.e., have large displacements of lifting surfaces predominantly in 
one direction; (2) have large motions of nonaerodynamic components with large 
mass, so that the generalized mass is high and the total percent of critical 
damping (aerodynamic plus structural) is low; (3) have large motions at critical 
vehicle locations, e.g. , crew compartments and payload area; and (4) are 
lightly damped aerodynamic ally due to proximity to a flutter condition. Of 
course, few if any modes will exhibit all these characteristics, and hence a 
good preliminary vibration and gust -response analysis will be required to 
delineate the "critical modes. Therefore, considerable knowledge about a 
vehicle's mass and structural characteristics must be available before a 
definitive selective judgement can be made along these lines. Generally 
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speaking, however, the implications of the above criteria are that a good 
simulation is required for: (l) the lifting surfaces; (2) the most flexible 

portions of the fuselages; ( 3 ) the interstage structure; (k) the crew- compartment , 
payload, and tank support structures; and ( 5 ) the engine and. sensor support 
structures. This simulation is required not only to provide the desired modal 
data, but also to provide data on internal structural loads due to inertia 
loadings associated with each mode. The simulation should emphasize those 
flexibilities that give rise to transverse rather than longitudinal deformations. 
Also, since liquid motions will give rise to large fractions of many of the 
generalized masses and internal loads, these must be represented carefully. 

6.3 Evaluation of Importance of Higher-Frequency Modes on Gust Response 

To obtain a preliminary quantitative evaluation of the number of modes for 
which data would be required to perform space-shuttle gust-response calculations 
of reasonable engineering accuracy, an external hydrogen tank orbiter heat 
sink booster, the Grumman/ Boeing H3T design in the ascent configuration 
corresponding to maximum dynamic pressure was analyzed using 22 of the first 
30 symmetric modes of the combined orbiter/booster. The 22 modes used were 
selected based on the criteria in the previous paragraph, and covered a 
frequency range from 1.5 to 8 Hz (see Table 6-1). The rigid body short-period 
mode was about 0.35 Hz. Frequency-response functions and power spectra were 
computed for the accelerations at 10 locations, 5 on the orbiter and 5 on the 
booster. Table 6-2 which summarizes the results, gives the order of importance 
of the. various significant modes according to their contributions to the mean 
square acceleration at each location. In attempting to correlate the peaks of. 
the power spectra in Appendix A with the mode numbers listed in Table 6-2, 
the shifts in the modal frequencies due to the aerodynamic forces must be 
considered. This was done with aid of modal response curves such as those 
shown in Figures A-5 through A- 8 . As might be expected, the first mode is by 
far the most important of the various structural modes considered. However, 
there are 11 different modes which have contributions greater than the rigid- 
body accelerations at least one location. - Furthermore,- the : 21st mode has the 
greatest contribution at two locations, and even the 29 th mode has a significant 
contribution at one location. 


In evaluating the significance of these findings from the viewpoint of 
dynamic modeling requirements, the limited scope of this preliminary study 
must be considered. Of the 22 modes chosen for inclusion in the calculations, 

13 were found to have significant contributions at least one of only ten 
locations. Thus, it is likely that there are other modes that would have 
significant contributions someplace on the vehicle. It appears, therefore, 
that the calculations performed thus far are approaching the definition of an 
upper limit, but that this limit cannot yet be stated with confidence. At 
this point, a best guess at a conclusion would be that a space shuttle structural 
dynamic model for gust-response purposes should probably be capable of rep- 
resenting modes up to about 8 Hz. A representation up to about 6.5 Hz is 
definitely required. 


TABLE 6-1 H3T COMBINED ORBITER/BOOSTER 


SYMMETRIC VIBRATION MODES 


MAX Q ALTHA 



25,240 

9,032 

569,811 

181,346 

51,781 


133,363 

18,355 

23,564 

5,215 

4,724 



37,682 

27,335 

56,611 

9,969 

201,599 


17,127 

■9,618 

41,576 

6,847 

2,890 


Mode 

Normalized 

to: 

Used in 
Gust- 
Response 
Analysis 

Orb Fin Z 

X 

Boost Wing Z 

X 

Orb Lo_ Z 

X 

Orb Lop Y 

X 

Orb Fin X 

X 

Orb LH Z 

X 

Orb Fin X 

X 

Orb: Fin X 

X 

Orb Wing Z 

X 

Boost Fin Z 

X 

Orb Fin X 

X 

Orb Wing Z 

X 

Orb LHg Y 

X 

Boost Fin X 

X 

Boost Wing X 


Orb Eng Z 


Boost Fin X 

X 

Boost Nose X 

X 

Boost Stab Z 


Boost Stab X 


Orb Nose Z 

X 

Orb Nose Z 

X 

Boost Wing Z 

X 

Orb Nose Z 

X 

Orb Lo^ Y 


Orb OMS X 
Orb Nose Z 
Boost Nose Z 
Orb Wing Z 
Orb Wing Z 

X 

X 
































TABLE 6-2 - RELATIVE MODAL CONTRIBUTIONS 
TO REPRESEHTATIVE LOCAL ACCELERATIONS . 
H3T COMBINED ORBITER/BOOSTER 


Location 

Order of Significant 
Modal Contributions 

Orbiter 



Crew Compartment, 

(Coord 57) 

21, lj RB *, 18 

Payload, 

( 60 ) 

21, 23, RB, 22, 10, 12, 29 

Oxygen Tank, 

(39) 

1, RB, 3 

Hydrogen Tank, 

(84) 

1, 21, 7, 22, RB, 23 

Wing Tip, 

(no) 

1, 9, 6 

Booster 



Crew Compartment, 

( 126 ) 

IB, RB, 7, 10, 3, 1 

Oxygen Tank, 

(150) 

RB, 1, 7, 21 

Hydrogen Tank, 

( 136 ) 

RB, 1, 7, 3, 21 

Wing Tip, 

( 173 ) 

- 2, 1, 22, 23 

Engine, 

( 196 ) 

1, 7, RB 


* RB :=‘ rigid-body short-period mode 
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7. ENGINE INDUCED RESPONSE 


7.1 Brief Description of Analysis Procedure 

Dynamic loads due to engine induced forces are generally determined by 
calculating the response of unrestrained vehicles represented by a selected set 
of free modes of vibration. Time integration is used for transient loading, 
and harmonic and random vibration analysis procedures are used for oscillatory 
loading conditions. The start-up forces for the shuttle will probably employ 
sequential ignition to reduce dynamic loads and can therefore be represented 
by a series of steps, each with a rise time of about 0.6 seconds for the liquid 
engines and 0.3 seconds for the solids. The time from 10 percent to full thrust 
is probably greater than .6 seconds. Engine shutdown forces will be more rapid, 
but the total sequence will probably occur over half a second. The SRM abort 
thrust termination could occur rapidly but the process could be controlled to 
occur over .2 sec. 

The random engine induced forces will cover a broad frequency range. 
Currently requirements limit the harmonic content of engine thrust to 450 
pounds for each liquid rocket engine in the frequency range from 1.5 to 25 Hz 
and permit slightly higher levels outside this band .up to 100 Hz. 

7.2 Structural Transfer Functions Required 

The transient and steady-state thrust irregularities occur primarily in 
the axial direction. Thus the strongest excitation arises for modes involving 
substantial axial engine motions. The largest overall loads derive from the 
lower frequency modes in which large sections of the vehicle move in phase, 
and which also involve significant engine axial motion. However, the 
acceleration response particularly on the orbiter tends to occur, at frequencies 
close to the engine support resonance. This effect is illustrated by the 
responses calculated on the Grumman GUI model of the Shuttle. Most of the 
higher frequency response occurred at 12.4 Hz, as shown on Table 7-1. In order 
to calculate responses and loads both, modes up to about 20 Hz would be required. 
If only loads were required, a lower cut-off frequency would.be acceptable. 
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TABLE 7-1 

ACCELERATION RESPONSE OF SHUTTLE CONFIGURATION GUI TO !t 50 LBS/ENGINE OSCILIATHC THRUST AT LIFTOFF 
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8. STAGING 

8.1 Staging Analysis 

Staging of the shuttle vehicle will normally occur at an altitude of 
approximately 130,000 feet and a velocity of about Mach 4. Under nominal 
staging conditions, the boosters are essentially empty while the liquid 
propellant tanks are about 2/3 full. In contrast, abort staging may occur at 
any intermediate altitude and associated velocity and propellant loading. The 
number of abort conditions is intrinsically large and will demand a larger 
share of the analysis effort. 

During nominal staging or in an abort staging, loads are input to the 
system from a variety of sources, all of which must be accounted for in the 
staging analysis. Initial conditions form a second input, to the staging 
analysis. The inputs and other factors must be considered in the analysis, 
including: 

• Stage separation systems which are kinematically constrained. These 
devices may. be pre-loaded and include spring elements or active devices 
such as pistons. 

• Thrust transients arising from shut down of booster engines. The 
booster engine thrust may not tail down uniformly which could give rise 
to differential thrust and require ME gimballing to counteract the 
motions. 

• Stored elastic energy in the orbiter and booster which is released upon 
staging. 

• Flight loads from aerodynamic and control effects. 

• Liquid propellant motions. 

The analytical structural dynamics model to determine staging response 
of the orbiter and booster will include the following elements: 

• Bending, longitudinal, and torsional deformations. 

• Slosh effects will be included implicitly in the modes if slosh 
amplitudes or masses are small. The large amplitude slosh motions which 
may arise during abort staging, such as propellant sidewall impact, can 
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be bounded using existing experimental data. A more precise deter- 
mination of this effect will depend- on data obtained from the structural 
dynamic model. 

• Detailed booster and orbiter thrust time -hi stories, including the 
thrust of each engine and the orientation of each gimballed thruster 
as a function of time. 

• Initial conditions such as angular rates, and velocities, which arise 
from pre-loads, detonation of pyrotechnics, and actuators. 

• Analytical models of the gasdynamic and aerodynamic effects. 

• Analytical model of stage separation system kinematics, stiffnesses, 
power sources, and release mechanisms. 

The above inputs will be combined in a numerical simulation of the stage 
separation process. The equations of motion will be synthesized from the 
equations of motion for the orbiter and booster written in terms of the flexible 
body modes and rigid body modes and the equations of motion of the link mechanism. 
The synthesization is achieved by invoking the boundary conditions between the 
stage separation system and the orbiter and between the booster and the stage 
separation system linkage. Because of the large angle motions of the booster 
and orbiter that may arise during some staging maneuvers, the large amplitude 
propellant motions, and the nonlinearity of the stage separation system, the 
equations of motion for the system will be nonlinear, necessitating a numerical 
integration of the equations of motion. 

Three major concepts for the stage separation system are presently being 
studied. These are a reverse linkage mechanism, a piston device, and 
combination, link-actuator system. Analyses have shown a major factor in 
establishing the forces and moments which arise during staging is the stiffness 
of the interconnection attachment points. Loads which arise during staging 
vary over orders of magnitude, depending on the stiffness of these points, 
and for some ranges of stiffness a "rocking" motion (pitch oscillation) can 
arise- during- staging.-- y - - ' ' - ' ~ “ “ ‘ ' . ' v 


The same considerations apply to abort staging. In addition, during 
abort staging (the worst conditions may exist at max qcr), large propellant 
loadings may be present. The fuel may accelerate rapidly and impact on the 
tops of the tanks, causing large upsetting moments and large bulkhead loads. 

At the present time, uncertainties in stage separation system design 
including stiffness, geometry (kinematics), release times, damping, and 
operation in conjunction with system control preclude a detailed modeling of 
the system. In any case, the properties of the vehicles at the points of 
attachment of the stage separation system linkage are the analytically uncertain 
quantities which are most readily established by means of the structural 
dynamic model. The actual stage separation system should be developed and 
tested by means of a separate structural scale model. 

The structural dynamic model should provide a portion of the inputs for 
the staging analysis: the modes and frequencies of the orbiter and booster 

separately, with coordinates defined at the stage separation system linkage 
attach points, and a means of characterizing the large amplitude slosh motion 
which are likely to arise during some abort staging modes. 

The experimentally established free -free modes and frequencies of the 
orbiter and booster will serve to validate the analytical procedures. These 
are used in conjunction with a separate dynamic model of the stage separation 
system and analytical models of loads and initial conditions to simulate 
staging dynamics. 

8.2 Structural Transfer Functions Required 

Depending on the design, the stage separation system will typically attach 
at 6 to 8 points.. Modal data must be obtained at these points on the booster 
and orbiter structural dynamic models. Accelerometers will be required to 
measure linear accelerations and angular velocity transducers will be used to 
measure modal rotations at these points. . 

Since staging will occur over a time period of the order of 1/2 sec., 
the low .frequency modes are most important to the staging analysis. A 
determination of all orbiter and booster, modes below 5 Hz will be adequate for 
staging purposes. 
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The stage separation system linkage itself need not be included in the 
structural dynamic model. However, since attachment point stiffness is of 
critical importance in determining staging dynamics, portions of the stage 
separation system which contribute significantly to the stiffness of the. 
linkage (or piston or actuator) attach points must be carefully, included in the 
structural dynamic model. 

A second area of major concern to the staging analysis is the modeling of 
the booster/orbiter thrust structures. Booster and orbiter angular and relative 
displacement rates which arise during the staging maneuver will be controlled 
to preclude structural damage to the booster and orbiter, .as .well as the 
constraints of crew safety and confort. Generally all booster and orbiter 
engines are thrusting during staging. Control is achieved by means of an 
appropriate transient variation in the thrust of the orbiter engines and by 
gimballing of the engines. Thus, flexibility of the thrust structures will 
strongly influence the angular and displacement rates each vehicle actually 
attains. Modal data must be defined for each engine (as also required on the 
basis of POGO analysis considerations) in. addition to interconnection attach- 
ment points. The discussion "Requirements for POGO Analysis," pertaining to 
thrust structure modeling are equally applicable for staging analysis purposes. 

Large amplitude sloshing and propellant impact play a major role in 
staging abort condition dynamics. In the case of nominal staging, the ullage 
of the External Tank could be about l/3. Significant propellants mass may 
therefore undergo large amplitude angular motions during the staging maneuver 
or impact tank bulkheads. Further study should be given to determine if 
separate tests will be required to dynamically characterize these motions, or' : 
it should be established by analysis that their occurrence can be precluded in 
every abort mode by means, of appropriate system control. , 

A study of the reverse link concept in an early Shuttle configuration has 

shown that during nominal staging the moment applied on the booster through 

6 

the stage separation system is approximately 4o x 10 ft-lbs, the moment produced 

6 

by plume impingement is -on the order of - 3 10 ft -lb, and the moment produced 

6 

by thrust is approximately 4o x 10° ft-lb. Significant variations in these 
number will occur before the stage separation system design is fixed and the 
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staging process is defined in detail, however, they provide order of magnitude 
indications of the magnitude of the forces required to study staging transients 
in the laboratory. 

The test configurations required for nominal staging are: 

Modal Survey 

• Free-free support of the booster with nominal staging equivalent 
propellant loading 

• Free -free support of the orbiter with full propellant loading 

The test configurations required for abort staging include: 

• Free -free support of the booster with propellant load corresponding 
to max qck' 

• Free-free support of the booster with intermediate propellant loadings, 
such as. 25 $, 50$ and 75$ 

Slosh Test 

A slosh test should be performed if the slosh response during abort 

staging can be shown by analysis to be large. These tests will consist 

of measuring vehicle and slosh response to simulate staging loads. 

If abort staging response tests are deemed necessary, it will have ! 
significant impact oh booster model test facility requirements. That is, this 
test may require special support fixtures, actuators for inputting simulated 
staging loads, and additional instrumentation to measure fluid and tank wall 
response and structural loads. 
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9. TEST REQUIREMENTS 





9. TEST REQUIREMENTS 

The test requirements for many of the analytical procedures are similar ; 
therefore this topic is discussed in this one section to limit repetition. 

9.1 Measurement Requirements and Instrumentation 

For most of the mode surveys the following measurement locations are, 
suggested. 


Orbiter 


Accel. 

External Tank 

Accel. 

Strain 

Pressure 

Wings 


20 

LOg Tank Wall 

12 

12 


Elevon 


5 

LO 2 Tank Dome 

10 

12 

1 

Fin 


8 

Interbank Skirt 
(SRM Attachment) 

6 



Fuselage 


8 

LH^ Tank Wall 




Crew Area 


3 

LHg Tank Frames 
(Orbiter Attachments) 

6 



Payload 


12 

Aft Skirt 
( SEM Attachment ) 

6 



Engine 


9 





Booster 


Accel. 





Nose 


4 





Fwd Mounting Skirt 

4 





Fuel Cylinder 


16 





Aft Mounting Skirt 
and Holddown 

6 






Instrumentation systems currently available, such as the automated mode 
survey data acquisition and reduction systems at Grumman ; and TEW are suitable, 
for most tests. The Grumman system for example has capability for up to 300 
measurements and can process accelerometer, strain gage, and force transducer 
signals. Automatic data plotting using computer assisted procedures are 
helpful in rapidly processing the data. Lightweight piezo electric accelerometers 
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suitable for model tests are available. Fluid reaction in LO^ tanks would be 
determined from the strain gages provided, sufficient sensitivity is available, 
alternatively, pressure transducers might be used. 

9.2 Facilities 

Facilities required for the model testing would include a low frequency 
suspension system to simulate free-free conditions, hoisting and handling 
facilities, tank fill and drain capability, and shaker control and support or 
suspension systems. For cantilevered tests, a base plate suitable for mounting 
hold-down latches is required. A shaker control system capable of operating 
8 simultaneous shakers should--be ; adequate. 

9.3 Additional Requirements 

Modifications or additions to the above test requirements are associated 
with data required for the following analyses. 

9.3.1 POGO - 

The significant data to be acquired is the response to axially applied 
excitation at the engine. Configurations to be tested are the combined vehicle 
at 3 different levels of propellant in the SRM booster, and the orbiter and 
external tanks at 3 different propellant . levels. Force inputs are required at 
the engine. Response of the fluid in the tank which can be used to determine 
tank bottom pressures should be measured. Data should be reduced to supply 
transfer functions between sinusoidal sweeps at the engine and response at the 
tank bottoms. 

9.3.2 Flight Controls 

Accelerometers and rate gyros should be . located in the vicinity of 
proposed Flight Control System sensor mounting locations as well as at the 
points of attachment of the control surfaces. Mode surveys should be conducted 
at 3 propellant weight levels for the combined vehicle and for the orbiter 
alone . 
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9.3.3 Ground. Wind Loads 


Cantilevered modes of the combined vehicle are required. Test config- 
urations probably should provide three values of hold-down latch stiffnesses. 

Four propellant configurations should be tested for each value of stiffness. 

9.3.4 Lift-Off 

Mode survey data for the combined vehicle on a soft suspension is the 
primary requirement. However, the response to transient inputs at the engine 
and hold-down points should also be measured. Techniques like the impulse 
transfer function tests would permit rapid and accurate calculation of the 
response to hold-down and engine transients which might be difficult to 
reproduce for a scaled model. 

9-3.5 Gust 

The most important configuration for which information is desired from 
a dynamic model for gust -response studies is the combined orbiter/booster in a 
vertical-flight position for a weight condition corresponding to the maximum- 
dynamic -pres sure portion of the ascent trajectory. A secondary objective 
would be to obtain similar information for about three other points during the 
portion of the ascent trajectory f ran lift-off to about 50,000 feet. In addition 
data would also be desirable for the orbiter independently for weight conditions 
corresponding to: (l) the maximum-dynamic-pressure portions of their re-entry 

trajectories; and (2) typical near-horizontal -flight configurations. 

9*3.6 Engine Induced Loads 

Test configurations should include the full combined vehicle, the maximum- 
dynamic -pres sure configuration, booster burn-out, and orbiter/external tank 
at 3 different propellant quantities. Single point excitation at the engines 
should be used to determine transfer functions for both steady and transient 
loads. Impulse transfer function techniques would be useful for determining 
responses to transient applied loads not readily reproducible in test. 


9.3.7 Staging 

The booster and orbiter would be tested separately in free suspensions. 
The configurations for the booster would include the appropriate propellant 
weights for lift-off, maximum dynamic pressure, and burnout. The orbiter/ 
external tank would only be tested with propellant weight at separation. 

Measurements of accelerations and rotations should be made at the inter- 
connection points in addition to the standard instrumentation. 

Slosh tests of the vehicle in maximum dynamic pressure configuration may 
be required in which simulated staging transient loads are applied and the 
resulting tank responses and slosh amplitudes are measured. 
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10. SUMMARY OF REQUIREMENTS STUDY 

10.1 Summary of Significant Structural Features for Each Dynamic Analysis 
Procedure 

The specific aspects of the shuttle structure which are most important in 
each of the. analyses procedures reviewed are summarized in this section. This 
information would be of direct use in designing an experimental model intended 
to simulate the structural dynamics of a specific prototype; however, it is 
also required to determine those features most significant and therefore likely 
to reveal problem areas. 

10.1.1 -P0G0 

• Tanks 

Sidewall and aft end dome bulging flexibility under fluid loading for 
L0^ tank. Flexibility of intertank skirt and LH^ tank to SRM and 
. orbiter interstage locations. 

• Engine 

Orbiter engine support structure, flexibility 

, • Feedlines ■ 

Not included 

• Orbiter Fuselage 

Flexibility from engine support structure to interstage 

• Booster 

Flexibility in lowest lateral mode 

• Frequency Range 

Modes involving motion of major parts of the launch configuration up 
to about 5 Hz are required. Modes involving motion of the orbiter 
engine thrust structure, feedline support points, and LOg tank wall 
and bottom dome are desired to about 25 Hz. 
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2 Control System 

• Engines 

Engine support structure flexibilities up to the first axial mode 

• Aerodynamic Control Surfaces 
Surface support flexibility 

• Reaction Control Thrust 

Flexibility of structure where thrusters are located 

• Sensor Locations 

Flexibility in areas where sensors might be mounted 

• Frequency Range 

Modes up to 7 Hz are desired for both the combined vehicle and orb iter 
alone . 

3 Ground Wind Loads 

• Launch Configurations 

Flexibility adequate to describe the first cantilever modes in bending 
. normal to the wing plane, lateral bending, and torsion (roll). 

• Aerodynamic Lifting Surfaces 

Wing flexibility adequate to describe the first symmetric mode 

• Frequency Range 

Modes up to 2 Hz are desired for the full vehicle restrained at the 
base by the hold-down latches. This same frequency range is also 
considered adequate for the empty vehicle. 

4 Lift-Off Loads 

• Fuselage 

Orbiter flexibility between the engine and the interstage points 

• Aerodynamic Lifting Surface 

Flexibility adequate to describe lowest two symmetric bending modes 

• Tanks 

Sidewall and aft end dome bulging flexibility . of LOg tank adequate to 
described lowest modes. Flexibility from LOg tank through LH^ tank to 
Orbiter and SRM interstage points 
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• Booster Thrust and Hold-Down Structure 
Flexibility equivalent to the first axial mode 

• Frequency Range 

Modes up to 4 Hz are desired for the full combined vehicle. 

5 Gust Loads 

• Wings 

Flexibility up, to the equivalent of the second bending and first 
torsion modes 

• Fuselage 

Orbiter flexibility equivalent to the first bending mode 

• Tanks 

Flexibility from L0 2 tank through Orbiter interstage points to wing 
root 

• Payload and Orbiter Crew Compartment 

The flexibility relative to the wing root 

• Booster 

Flexibility equivalent to 1st axial mode 

• Frequency Range 

Modes up to about 8 Hz are desired for the combined vehicle in maximum 
qo? weight configuration. 

6 Engine Induced Response 

The same characteristics as required for POGO analyses would be suitable 
except that the frequency range of interest would extend somewhat higher. 

7 Staging 

• Engines 

Engine support structure flexibilities up to the first axial mode for 
both the orbiter and the booster 

• Interstage Connection Points - 

Flexibilities from these points to major inertia reaction items such as 
the L0 2 tanks 
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Frequency Range ■■ 

Modes up to 5 Hz for each vehicle component separately -for a full 
orbiter and for an empty booster. 



11. MODAL COUPLING APPLIED TO THE NASA LANGLEY 1/15 SCALE 
STRUCTURAL DYNAMICS MODEL TEST DATA 



11. MODAL COUPLING APPLIED TO l/l5 SCALE MODEL TEST DATA 

In accordance with the requirements of NASl-10635-^ Modification 1, 
analytical coupling of the orbiter and booster dynamic model, fuselage modes 
has been briefly investigated. This work was undertaken to evaluate a specific 
coupling procedure using test data from the 1/15 scale model described in 
References 11-1 and 11 - 2 . 

11.1 Description of Problem 

An additional requirement to this task added by Modification No. 1 on 
9/1/71 stated that "In order to more fully determine the requirements for the 
model configuration, additional study shall be conducted to further evaluate 
modal component synthesis techniques needed to attain an adequate understanding 
of the dynamic behavior of the involved structure. " The objective was to 
evaluate an existing analytical approach and to study means for computing 
systems damping from component tests. Therefore the feasibility of analytically 
combining the experimentally determined component modes of a model of an early 
shuttle concept was investigated. 

The model, described in Reference 11-1 and 11-2, consisted of tubular 
beams representing an orbiter and a booster in pitch plane motion. These beams 
are coupled by two interconnecting spring assemblies each of which had a 
stiffness of 1,000 lbs/in.- The measured data and analysis only considered 
lateral motions. 

11.2 Description of Coupling Procedure 

The. modal coupling procedure used assumed that the orbiter and booster 
fuselage models could be coupled by massless springs which add to the generalized 
stiffness of the combined beams. The potential energy of the coupled system 
would then have a contribution from the orbiter, described in terms of similar 
modes, and a contribution from the coupling springs between the two vehicle. 

It was assumed that the kinetic energy of the coupled system could be described 
terms of uncoupled orbiter and booster modes. The equations of motion of the 
combined system, in terms of model coordinates as determined from the energy 
expressions, were solved for the eigenvalues and eigenvectors, and the resulting 
vectors were translated back into the original coordinates. The following 
is a summary of the procedure. 
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For the specific problem considered here, coupling 2 flexibile and 2 rigid 
orbiter modes to 5 flexible and 2 rigid booster modes, the size of the matrices 
involved are shown. Also, for equal rate coupling springs, the form of (K c ) is: 
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Similarly, the term required for the equations of motion from the kinetic 
energy becomes: 
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11.3 Verification of Coupling Procedure 


To verify the procedure, the simple two-beam, three -spring, arrangement 
shown below was examined. Each beam was allowed 14 degrees of freedom and had 
a fundamental (free-free) frequency of about 86 Hz. Two sets of equal rate 
springs were studied (10^ and 10^ lbs/in.). 
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A review of the calculations indicated that when l4 modes were used to 
describe each beam (two rigid body and 12 flexible) all 28 coupled-system 
frequencies agreed with those calculated directly for the entire system. The 
frequency comparisons are included as Tables 11-1 and 11-2. The effect of 
using fewer modes in the coupling process was also explored and it was determined 
that reproducibility of frequencies appeared to be better with a higher number 
of modes. For example, coupling the 10 lowest frequency modes of each beam 
gave correct total system frequencies up to the 20th mode; whereas , coupling 
3 modes of each beam gave only 3 correct total system frequencies. This 
result was influenced by the value of the coupling spring and generalizations 
do not appear warranted. Mode shapes are also compared in Table 11-3. 

The agreement shown in the above problem provided the desired check of 
the analytical method. 1 * 

11.4 Processing of Data from Langley 1/15 Scale Model 

The separate modes of the orbiter and booster fuselage models obtained by 
test in October , 1971 were used in the coupling analyses. The 2 elastic and 
2 rigid body orbiter modes used are shown in Table 11-4. Only motion in the 
pitch plane was considered. It was most convenient to keep the modal displace- 
ments at the locations of the concentrated weights separate from those for the 
distributed weights, therefore the first 7 coordinates are for distributed 
weight locations and the succeeding 16 are for the concentrated weights. The 
distance alongthe beam, from the orbiter nose, for each coordinate may be 
obtained by subtracting 8 in. (the reference distance) from the value for 
mode 4, the rigid body pitch mode. The weight associated with each (coordinate 
is shown in the column matrix on the right under WEIGHT 0. In order to check 
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2 

the orthogonality of these modes, the matrix of terms M. /M M was formed 

ip 1 J H J J 

from the generalized mass matrix M = (cp) (m) (cp). In evaluating experimentally 
measured modes, the orthogonality between them is generally considered adequate 
if the off-diagonal terms in this matrix are less than 0.1. This is the case 
for the orb iter elastic modes, the first two rows and columns in matrix ORTHO 
shown on the bottom of Table 11-4. The last two rows and columns show the 
coupling between rigid body translation and pitch modes which is high as 
anticipated because the rigid body pitch is taken about a reference point 8" 
forward of the nose. 

The 5 elastic and 2 rigid body booster modes used are shown in Table 11-5. 
The first 13 coordinates represent the distributed weight location, while the 
remaining 18 represent the booster concentrated weights. The weight associated 
with each coordinate is shown in the last column on the right under WEIGHT 8. 

The distance of each coordinate from the booster nose reference point may be 
obtained by subtracting 15.7 from the values in the 7th mode (rigid body pitch). 
The orthogonality check shown at the bottom indicates that only orthogonality 
between the 4th and 5th elastic modes is not within the desirable limit. 

The generalized mass matrix was formed from the measured flexible and the 
computed rigid body mode shapes. All off diagonal terms (except for the 
coupling between rigid body translation and pitch) were eliminated. The form 
of the matrix used is shown in Table 11-6. 

The generalized stiffness matrix was formed as indicated previously. The 
off-diagonal terms were eliminated from the separate body generalized mass 
matrix before multiplying by the square of the resonant frequencies. The 
coupling generalized spring matrix which was a full 11 by 11, was then added. 

The resulting 11 x 11 eignevalue problem was then solved and the resulting" 
eigenvectors were translated into the original booster and orbiter coordinate 
system, and then rearranged for easier plotting. Tabulated values of the 
modes and frequencies of the coupled system are shown on Table 11-7. The 
frequencies are shown at the top of each vector column. The location of each 
coordinate -may be -obtained from the 0.19 Hz mode” which turns out to be the 
rigid body pitch about some point just forward of the orbiter. The non -zero 
frequency of this rigid body mode results from the numerics of the eigenvalue 
calculation. The components of this modal vector represent the relative 
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distance of each coordinate, and multiplying them by 150 in. gives the absolute 
distance. For example, coordinate 15 is .517 of the distance from the reference 
point to the back end of the beam representing the booster, or 150 (. 517 ) = 78 ". 
The first 31 coordinates represent the booster, the next 23 represent the 
orbiter and the last 4 represent the location of the coupling springs on the 
orbiter and booster. 

11.5 Review of Results 

Five of the flexible modes for which comparisons between analytical 
coupling and direct measurement could be made were plotted as shown in 
Figures 11-2 through 11-6. One mode measured experimentally, the 4th flexible 
at 48.3 Hz did not appear at all in the analytically coupled modes but is shown 
in Figure 11-7. This mode shape resembles that for the second flexible mode 
with considerably more orbiter bending. An orthogonality check of the measured 
coupled system modes, as shown in Table 8, indicated that the 2nd and 4th 
flexible modes were highly coupled (note that the 2, 4 term in the matrix is 
O. 85 ). Both modes are apparently orthogonal with all the others measured. 

Modes above IO 3.5 Hz were not plotted because a review of the data in 
Table 11-7 indicated that the 4 remaining higher frequency modes were similar 
to the input modes for the individual uncoupled model in both frequency and 
shape. In fact, the sixth measured flexible coupled mode shown in Figure 6 
actually appears to be the first orbiter bending with little motion of the 
booster and the measured frequency of 101. 9 Hz is close to the 98 . 8 . Hz measured 
for the orbiter model alone. This effective uncoupling of higher frequency 
modes is due to the relatively low stiffness of the connecting springs. A 
similar effect was noted in the analytical work described in reference 11-2. 

The coupled pitch mode at 104.9 Hz shown in Figure 22(i) of this reference is 
similar to the first mode of the orbiter fuselage at 103.5 Hz shown in Figure 13 
(a) of the same reference. Also, the next higher pitch mode at 118.7 Hz shown 
in Figure 22(1) is similar to the second mode of. the booster fuselage at 117.7 Hz 
shown in Figure 12(b). 
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11.6 Discussion and Recommendation 


The coupling procedure employed in the investigation demonstrated the 
feasibility of predicting coupled system modes using vibration measurements 
made during tests on the separate vehicles. The inability of this procedure 
to predict the fourth mode (48.3 Hz), however, demonstrates the difficulties in 
making assumptions as to the type of restraint at interfaces between vehicles. 

In subsequent work reported by Dr. R. W. Fralich (Reference 11-3), the 
assumption that the interconnecting springs were pinned to the tubes was 
modified to permit a finite rotational restraint, and the assumption of lateral 
motion only was modified to permit axial motion as- well. With, these changes, 
the calculated and measured frequencies for the fourth as well as the other 
modes were very close. 

The case studied here involved physical models of relatively simple form, 
i.e., two stiff beams represented by these free-free modes coupled by relatively 
soft discrete springs. The evaluation of other coupling techniques which 
might be used on the shuttle should consider more representative model con- 
figuration and a wider variety of testing and coupling combinations. The 
existing model, if it is to be used further, could be modified to account for 
more rigid interstage connections (including redundant ties) and for motions 
other than in the pitch plane. Close coordination between analyses and test 
should be maintained to insure the careful investigation of . discrepancies. 
Analytical coupling procedures which accommodate separate vehicle tests under 
other than free -free support conditions should also be investigated. 
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SAMPLE BEAM PROBLEM 

COMPARISON OF MODE SHAPES FOR ANALYTICAL COUPLING 
K = lb S /fh 


Mode Freq Mode Shape Calculated with Number of Modes Shown 

No. (KPS) Agrees with that From Full 28 Coordinate Calculations 
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K = 10 


n 





3 

542 

0 

0 
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NOTE: 0 indicates aggreement 

X, indicates non-agreement - - . 


TABLE M-3 
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12. DESIGN OF 1/10 SCALE PRELIMINARY MODEL OF SERIES - BURN SHUTTLE CONFIGURATION 

12.1 External Configuration and General Arrangement 

After the series-burn configuration had been selected as the baseline, 
Modification 3 to this Task Order was issued to change the remainder of the 
effort at GAC to study the design of an appropriate structural dynamics model. 

In particular, fabrication costs and time were to be determined and the details 
of the propellant storage system were to be studied. 

The external configuration of the model is shown schematically on Figure 
12-1, The design details studied for the various parts are shown on the drawing 
numbers called out near each part. These drawings have been previously 
submitted for NASA review. Reproducible copies of these drawings have been 
sent to the NASA Langley Research Center Technical Monitor for this task. 

A general description of the model construction for the Orbiter, Booster, 
and HO tank assembly is given in the Stress Analysis Report in Appendix B. 

12.2 Details of the Propellant Storage System 

The L0 2 tanks were of welded aluminum monocoque design. The HO tank was 
.032" thick within 1/2" of the welds and chem milled to .016" over the remainder 
of the surface. The dome was elliptical of .016" thick aluminum. The Y ring 
was designed as a stepped element with provisions for welding the dame and the 
cylindrical portion of the tank. The booster L0 2 tank was heavy walled (.1" 
dome, .167" cylindrical section) aluminum. Simple flat cover plates were used 
in place of upper domes on both tanks. 

The LH 2 tank was a thin walled riveted tank, with a .008" thick aft dome 
and a cylindrical section which was .002" thick for the 120 degree segment 
nearest the orbiter and .016" elsewhere. In order to stabilize this structure 
against 1.5 g limit handling loads, rings were, added at about 2 1/2" pitch. 

The Y ring was machined in a stepped section. 

The RP-1 tank at the base of the booster consisted of an aft elliptical 
end dome of .113" thickness, a machined Y ring, a cylindrical section of .139" 
thickness, a forward machined Y ring, and a forward dome with access through an 
elliptical cover plate. This tank was also seam welded of aluminum. 
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12.3 Interbank Structure 

The skirt between the LOg and LH 2 HO tanks was .15" thick skin stiffened 
by angle rings. At the center of the sjcirt two back-to-back channels formed a 
heavy frame with the fittings for the orbiter attachment fastened between them. 

A similar arrangement was used inside the LHg tank as a support for the aft 
fitting from the orbiter.. The members between orbiter and HO tank were simple 
tubular elements. The skirt between the HO tank and the booster was of similar 
ring stiffened .015" thick skin. The aft thrust structure on the booster was a 
conical element with 12 radial bulkheads meeting in a hexagonal shaped cylinder 
in the center. ' - 

12.4 Orbiter 

The orbiter mid fuselage was basically series of U shaped frames with 
corner longerons along the top and bottom covered by aluminum skin. The skins 
were stiffened by intermediate frames every 2 " to prevent buckling. The orbiter 
engine support structure consisted of full bulkheads connected by longerons. 

The wings consisted of four beams carried through the fuselage covered by flat 
skins. 

12.5 Budgetary Cost 

A budgetary cost estimate was performed by a member of the manufacturing 
group normally making shop estimates at Grumman. Preliminary drawings AD6o4 
through 611 defining the HO tank and booster were used. The 1970 rates were 
applied including G&A costs but not including profit. Results were $55>380 
for the HO tank and $64,200 for the booster. The cost of the orbiter was not 
estimated. Fabrication- time was estimated as 4 months from go-ahead. 

12.6 Utilization of Designs 

The design work on the orbiter and the external tanks completed in this 
contract was used subsequently as the basis for the, design of a 1/8 scale 
structural dynamics model of the parallel burn shuttle using solid rocket 
boosters. 
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APPENDIX A 

A. 

GUST RESPONSE ANALYSIS 


A.l 

Description 



In order to provide representative quantitative data on the response of 
a shuttle type of vehicle to gusts, an extended and improved version of the 
aircraft program described in Reference (A-l) was utilized. In this program, 
vehicle motion is described by a linear, multi -degree -of-freedom system of 
forced response equations in terms of rigid and flexible normal modes. The 
effect of control laws in deflecting the engine due to rotation and rotation 
rate is included in the two rigid-mode equations. Additional differential 
equations, derived from an indicial formulation of the unsteady aerodynamics, 
are used to define the generalized aerodynamic forces due to motion and gust. 
The program can consider two symmetric rigid-body modes (vertical translation 
and pitch) and up to 8 symmetric vibration modes of the wings and fuselages. 

In order to accommodate the 22 flexible modes which appeared significant, the 
program was run three times, with the first flexible mode appearing in all 
three runs. Although this is not rigorously correct due to the elimination of 
the aerodynamic coupling between many of the modes, it is considered acceptable 
for this application. The modes used in each case are as follows: 


Run 1 

(15 July 1971) 

Run 2 

• (23 July 1971) 

Run 3 

(29 July 1971) 

Mode 

Freq. 

Mode 

Freq. 

Mode 

Freq. 

No. 

(Rad/Sec) 

No. 

(Rad/Sec) 

. No. 

(Rad/Sec) 

1 

9.23 

1 

9.23 

1 

9.23 

2 

9.84 

6 . 

15.77 

3 

13.54 

9 

18.78 

7 

17.14 

4 

14.52 

11 

20.25 

10 

19.26 

5 

. 14.89 

. 12 

20.59 

21 

34.99 

8 

18.25 

17 

27.46 

23 

39.12 ' 

13 

22.28 

22 

38.52 

29 

48.78 

14 

23.81 

24 

4o.51 

30 

50.23 

18 • 

‘ 

29.07 
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In formulating the generalized forces, strip theory, appropriately 
modified to yield experimental complete -vehicle steady-state stability deriva- 
tives, was used. The orbiter wing and booster wing were each represented by 
6 panels per side, with the aerodynamic loading on the remainder of the two 
vehicles being assumed to act on center section panels of each of the two wings 
Values of the coefficients used are shown on Table A-l. 

The shuttle configuration used in the analysis was an external hydrogen 
tank orbiter on a heat sink booster designated H3T, as described in Reference 
(5-2). The vehicle was assumed in steady flight at Mach 1„05 at 25,000 feet. 

The lowest 30 symmetric modes for the ascent configuration corresponding 
to maximum qa were calculated, and the 22 of these which came closest to 
satisfying the criteria of Section 6.2 were used in the gust response analysis. 
The general characteristics of these modes are described in Table A-2. 
Structural damping was assumed as 2 percent in each mode, after initial analyse 
showed that the aerodynamic damping in some modes could be significantly lower. 

The turbulence was represented by the boost spectrum shown in NASA 
TM X- 53973, Section 1.3.2., 

A. 2 Results 

The analysis proceeds by calculating the response to a 1 ft/sec sinusoidal 
gust velocity applied at various frequencies in the range of interest. The 
transfer function obtained is then used in calculating the response to the 
spectrum representing the turbulence. For each mode and coordinate therefore 
both a frequency response to a sinusoidal input and a power spectral density 
response is determined. 

The generalized force frequency response illustrates the effects of 
unsteady aerodynamics including wing sweep and transport time delays between 
the two surfaces. The power spectrum of the generalized force exhibits a 
single low frequency peak at about 11/2 radians per second, reflecting the 
nature of the turbulence which peaks at about 2 radians per second and then 
drops off at about 20 db per decade. These characteristics are shown typically 
on Figures A-l to A-4. 


A-2 



The modal accelerations in most flexible modes (Figures A-6 and A-8) show 
responses typical of uncoupled single degree of freedom systems; however, some 
coupling between modes is evident in the rigid body translation plots in 
Figure A-5 and in one higher flexible mode shown in Figure A-7. The integrated 
response to the applied turbulent spectrum yields the rms values shown in 
Table A-3. It may be noted that both the force and acceleration values remain 
significantly high even among the higher modes. 

The important feature in determining the significance of the higher modes 
is the contribution of such modes to the rms response at various locations. 
Responses were calculated for 5 locations in the orbiter and 5 in the booster. 
These are summarized in Table A-4. An estimate of the peak acceleration may be 
obtained by multiplying the rms values shown by 3. The largest acceleration in 
the orbiter coordinates other than the wing tip, occurs at the LHg tank and is 
equivalent to 0.45g. The crew compartment response is somewhat lower (0.4g)./ 
Responses in the booster are all lower. The contribution of each mode to this 
rms response may be determined by comparing the energy for the various modes 
in the power spectrum response at each coordinate. As an example Figures A-9 
through A-ll are the accelerations at coordinate 6o, the orbiter payload. 

Figure A-9 includes the 8 flexible modes used for Run 1 on 15 July, whereas the 
next two figures include the flexible modes used in Runs 2 and 3 on 23 and 29 
July, respectively. The largest portion of the energy under the power spectrum 
curve in Figure A-9 is at 2 1/2 radians/sec. , the rigid body (RB) mode. The 
energy near 38 radians/sec, the 22nd mode, is almost as high, and the mode at 
20 radians/second, the 12th mode, also contributes significantly. Figure A-10 
indicates that the flexible modes selected for the 23 July run were more 
significant, in particular, both the 21st and 23rd modes near 35 radians/sec 
encompassed much more energy than the rigid body mode, and the 10th and 17 
radians/sec and the 29th at 49 radians/sec are both marginally significant. 

The power spectrum curve on Figure A-ll indicates that none of the flexible 
body modes used in the 29 July run were as significant as the rigid body (RB) 
mode. 



Reviewing the power spectrum curves at other coordinates, it can be 
noted that the first flexible mode at 9 radians/sec is most significant for 
the remainder of the orbiter locations, except for the crew compartment (57) 
as shown in Figure A-l6, where the 21st mode contains more energy. Coordinate 
126, the booster crew compartment, responded in the 18 mode at 29 radians/second 
to a greater extent than in the RB mode. 

Response at each coordinate were reviewed in a similar manner and the 
results summarized as shown on Table 6-2 ou page 6-5. 


Table A-l 

Distribution of Lift Curve Slope Coefficients 


Panel 

No. 

Orbiter 

- ■ Booster 

Width 

(in.) 

• C L 

a 

Width 

(in.) 

C L 

O' 

.1 

185 

1.83 ■ 

190 

4.09 

.2 

80 

2.43 

200 

4.35 

.. 3 

90 

2.61 

200 

4.51 

4 

8° 

2.80 

180 

4.83 

5 

8° 

3.00 

200 

5.10 

6 

130 

2.58 

200 

3.43 


NOTE: 1. The panel 1 values represent the orbiter fuselage and 

the booster fuselage and horizontal tail. 

2. The integrated lift and moment curve slopes of the 

combined orbiter and booster are C_ = .6.55 and 

L_ 

or. 

= 1,1+9 compared with 6.65 and 1.77 estimated from 
<y • - 

wind tunnel model data on somewhat similar models. 







Summary of H3T Max Go 1 Mode Patterns 
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4.13 

4.37 1st Bending 1st F & A, : 2nd Bending, 2nd Torsion, 1st F & A, 

2nd Bending 1st F & A 1 1st Bending 2nd Bending 

4.63 1st Bending 2nd Bending 2nd Bending 2nd Torsion 1st F & A, 

2nd Bending 
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Table A- 3 


Mode 

Generalized 

Generalized 


Force 

Acceleration 


RMS lbs. 

RMS ft/sec 2 


Rigid Translation 

69>30 

1.09 

Rigid Rotation 

9,186 

• 52 


1 

(1.47Hz) 

9,288 

24.23 

2 

(1.57 Hz) 

11, 370 

21.96 

3 

(2.16 Hz) 

12,080 

1.15 

4 

(2.31 Hz) 

H9 

.12 

5 

(2.37 Hz) 

1,371 

1.49 

6 

(2.51 Hz) 

8,584 

3.39 

7 

(2.73 Hz) 

5,335 

6.75 

8 

(2.90 Hz) 

517 

2.06 

9 

(2.99 Hz) 

3,003 

17.13 

10 

(3.06 Hz) 

1,375 

8.34 

11 

(3.22 Hz) 

4,674 

8.42 

12 

(3.28 Hz) 

1,350 

2.43 

13 

(3.54 Hz) 

780 

• 51 

14 

(3-79 Hz) 

248 

1.47 

IT 

(4.37 Hz) 

5,153 

1.36 

18 

(4.63 Hz) 

4,460 

3.87 

21 

(5-57 Hz) 

8,661 

6.97 

22 

(6.13 Hz) 

10,680 

8.63 

23 

(6.22 Hz) 

11,430 

:4.58 

24 

(6.45 Hz) 

2,350 

3.05 

29 

(7.76 Hz) 

734 

2.95 

30 

(7-99 Hz) 

4,013 

12.18 
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STRUCTURAL ANALYSIS 

Structural analysis for the l/lO scale model design was limited to ■ 
checking the critical areas for buckling under a 1.5 g limit axial load with 
the vehicle in an upright launch attitude. A factor of safety of 1.5 to 
ultimate load was used. The stress analysis effort consisted of determining 
the approximate number and location of stiffening frames or stringers 
required. No additional stress calculations were considered necessary 
because the basic dimensions of the model were scaled down from the prototype 
which had been designed for critical flight loads and would therefore be 
adequate for any of the _ relatively ; low vibration loads induced during mode 
surveys. In view of the preliminary nature of the model design, no more 
detailed stress analyses appeared justifiable. 

An overall static equilibruim load distribution of the mated configura- 
tion under the ultimate 2.25 g load condition is shown below. 
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Model Configuration 

ms 1-10635-^ 

STRESS ANALYSIS 

H-0 External Tanks 

The oxidizer tank utilizes monocoque type construction. The area aft 
of x station 50.0 is critical for compression buckling and requires intermediate 
frames spaced between 2 £ and 3 5/8 in. This area includes the fwd. intertank 
skirt, LHg tank, and aft interstage skirt. Major bulkheads are located at the 
fwd. interstage location (x57*5), the aft interstage location (xl30.4) and the 
lag redistribution frame (x!12 .7 ) . Basic skin is aluminum alloy, .015" thick; 
it is locally thicker aft of the fwd. interstage location in the upper quadrant 
for increased buckling capability due to compression bending. 

Booster 

The booster I0 2 and RP-1 tanks are of monocoque structure construction. 
Thickness of the I£>2 tank is .167 in., and the RP-1 tank thickness is .139 in. 
The intertank structure aft of (x270.0) and fwd. of (x302.0) requires no 
intermediate frames and has an effective thickness of .0^1 in. The skirt 
aft of (x329.7) which houses the engine support beams is effectively .025" 
thick, but requires intermediate frames and stiffeners to prevent buckling. 

There are Seven main engines, one in the center and six outer engines which 
are supported by twelve support beams that run between a central six-sided 
cannister that houses the center engine and the outer shell. 
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STRESS ANALYSIS 


OKBITER / 

The orbiter fuselage between x57.0 and X130.40 basically consists 
of thin gage aluminum sheet with corner longerons made up of angle sections 
upper and lower. The skins are stiffened by intermediate frames at a pitch 
of approximately 2" with major bulkheads at the forward interstage fitting 
(x57-50) and the aft interstage location (xl30.4o). The forward portion of 
the fuselage is to be mass simulated as is the wing structure. The aft end 
of the fuselage model shall be simulated in such a manner as to effectively 
reproduce the true mass and stiffness distribution of the engines, support 
structure and tail structure. 
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